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Abstract 

We investigate in some detail the structure of the electromagnetic 
current density for the pseudo-relativistic massless spinor effective 
model for graphene. It is shown that the pseudo-relativistic massless 
Dirac field theory in 2+1 space-time dimensions and in the presence of 
a constant homogeneous electric field actually leads to the measured 
current density and to the minimum quantum conductivity. 

PACS numbers: ll.lO.Wx, 02.30.Sa, 73.43.-f 

1 Introduction 

One among the many far intringuing features of graphene is that [1] its 
conductivity at zero magnetic field doesn't vanish close to the charge carriers 
neutrality point of concentration n < 10^ cm~^ . Actually, on the one hand, 
the measured values for the lowest conductivity (Xmin near the neutrality 
point are very close to the conductivity quantum e^/h per carrier type |2j. 
This zero-field quantum conductivity does not apparently depend on the 
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chemical potential and keeps its single carrier value down to liquid helium 
temperatures [IJ. On the other hand, it turns out that the insofar known 
theoretical predictions fail [3J in reproducing (Jmin, a disagreement known as 
the mystery of the missing pie, but for a very recent calculation [4j. There, 
however, the deduction for the actual value of amin is certainly elegant and 
correct but admittedly rather formal. In particular, the direction in the 
graphene's sample plane of the minimum quantum current density is not at 
all clear from the above mentioned theoretical derivation. 

As a further important ambipolar electric field effect, which was not explained 
in |4j, it turns out that the single layer graphene resistivity p decreases rapidly 
to zero with adding charge carriers, showing their very high-mobility pQ. 
It is the aim of the present short note to fulfill this gap, by developing a 
new derivation for the canonical quantum current density, within the well 
established effective model [H [9] of a massless Dirac field theory in 2+1 
space-time dimensions and in the presence of constant and homogeneous 
electromagnetic classical background fields. In fact, the latter field theoretic 
model appears to be exactly solvable, in such a manner that all the gauge 
invariant physical observables, like e.g. the mimimum quantum conductivity, 
can be explicitly calculated and eventually compared with experimental data. 

In ref. [1] it has been shown that the Euclidean effective action, or the grand 
potential in statistical mechanics language, in the presence of a background 
uniform electrostatic field turns out to be independent from the (complex) 
chemical potential and from the temperature. Conversely, in the presence 
of a uniform magnetic field, the dependence of the effective action upon the 
chemical potential is crucial in order to reproduce the measured values for 
the quantum Hall conductivity [S] . It follows thereby that in the present note 
we can restrict ourselves, without any loss of generality, to zero temperature 
and zero chemical potential planar quantum electrodynamics in a constant 
homogeneous electric field. Hence, according to the leading order expansion 
around the Dirac points at the corner of the Brillouin zone in graphene [5], the 
classical pseudo-relativistic lagrangian in the Nikishov [TD] temporal gauge 
= (0, -Ect,0) reads 

C = ^{x)(^h^^ + -A^^^{x) ij = ij^ai (1) 

V c J vf 

Here ( — e) is the electron charge, c the light velocity in vacuum, vp ^ 1/300 
is the Fermi velocity of the massless Dirac quasiparticles in graphene, are 
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the rescaled Dirac matrices, i.e., 



7° = C7° ^' = vfY (^=1,2) 

while Fqi = Ex = E > denotes the constant and homogeneous electric field 
pointing towards the positive OX axis. We remark that in the planar QED 
on the 2+1 dimensional Minkowski space-time the spinor field has canonical 
dimensions of the inverse of a length, i.e. of a wave number. 



2 Uniform electric field 

Consider now a graphene sample in the presence of a constant electrostatic 
field pointing towards the positive OX axis, i.e., 

Foi = = = {E>o) 

In this section, we will choose a gauge which leads to non-stationary sets 
of solutions. In particular, to solve the Dirac equation in the present 2+1 
dimensional massless case, it is convenient to employ the representation for 
the gamma matrices 

7° = f^i 7^ = ^cr2 7^ = icTs (2) 

After setting = {ct,x,y) we get the massless Dirac operator in 

the Nikishov temporal gauge ^4^ = (0, Ect, 0) , i.e., 

iU 

\p— — dt(Ti + {ih dx + eEt) ia2 — hdyUz (3) 
Vf 

It follows that in the Nikishov temporal gauge the 1-particle hamiltonian 

Ht = fiVF [ (73(i dx - eEt/fi) - i 9^, ^2 ] (4) 

is explicitly time dependent, which means that the energy does not lead to 
a good quantum number for the specification of the Fock space states of the 
quantized massless Dirac spinor field. Thus, at variance with the magnetic 
field case, in the present context the Dirac hamiltonian does not play here 
any key role. Notice that The classical Action for the massless Dirac field in 
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the 2+1 dimensional Minkowski space-time and in the presence of a uniform 
background electric field takes the form 



/oo /• _ 

dt J dr *(t,r) 



— dfCTi + {ih dx + eEt) ia2 — hdya^ 



whence it is apparent that the canonical physical dimensions of the spinor 
field are equal to the inverse of a length. In order to obtain the solutions in 
this gauge, it is convenient to introduce the spatial Fourier transforms 



1 rOO POO ^ 

= — / dk dp e^P'^+i^y ^{t,p,k) 

ZTT J-oo J-QO 

with p = {px-iPy) = (p, k) as well as the quantum electric length 



= 

and the dimensionless quantities 

/ eEt 
C = [P- 



'eE 



n 



A = k''£l 



so that the Dirac operator can be recast in the suitable form 



— ik£E 



ikis 



If we set 



^{t,p,k) 



ip{t,p, k) 
x{t,P, k) 



we obtain the coupled differential equations 



{-id^ + ^)ip + ik£EX^'^ 



Then, we can write 



X = ( i/klE) {-id^ + ^ 
+ ^2 + A + i) = 



MO) 



(5) 



(6) 



(7) 



(8) 



(9) 



(10) 



k = 0) 



(11) 



Thus, in the presence of a homogeneous electric static field, one actually finds 
the following complete and orthonormal sets of time dependent solutions of 
the Dirac equation : namely, for A; 7^ we have either the outgoing normal 
modes 

1 



exp |i p ■ r — I vrAj 



(12) 



^ ' l{l+^)^/\D,x/2-li-Z-] 

[ Dix/2{-Z_) 

which correspond to the normal modes of an outgoing quasiparticle (massless 
electron or negatively charged neutrino) with momentum p and charge — e , 
while 

1 



27Tif 



exp |z p ■ r — 1 7rA| 



X 



iX/2 



(-^ + ) 



1(1 



(13) 



do correspond to the normal modes of an outgoing antiquasiparticle (hole) 
with momentum p and charge e , where we have set 

z± = (l±Oe 

Furthermore the incoming normal modes will be given by 

1 



(i,r) 



271 if 



exp |z p ■ r — I vrAj 



X 



-l{l+^)^/\D,x/2-l{z-) ^ 
Di\/2 {z-) 



(14) 



which correspond to the normal modes of an incoming quasiparticle (massless 
electron or negatively charged neutrino) with momentum p and charge — e , 
while 

1 



(i,r) 



exp I? p ■ r — I vrAj 



X 



(^- 



i\/2 



\/2-l[.Z + ] 



(15) 



^Thc complete and orthonormal sets of solutions listed below do coincide with those 
ones obtained in up to the overall dimensional factor ^ . 
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do correspond to the normal modes of an incoming antiquasiparticle (hole) 
with momentum p and charge e . 

Finally, for py = k — we come to the so called longitudinal normal 
modes or zero modes 



Up± {t,x) 



1 



27r£ 



E 



exp!^ipx+ ^C^(^)} 5 = Up{t,x) (16) 



Vp± {t,x) 



/ 1 \ 



1 




= Vp{t,x) (17) 



The above sets of time dependent massless Dirac spinors have canonical 
physical dimensions of a wave number, as it does, and are normalized in 
order to satisfy the following orthonormality relations, viz.. 



j dr u}^± {t, r) Up± {t, r) = i^'^ 6{p - q) 
j dr vl± (t, r) Vp± (t, r) = S{p - q) 
J drul^±{t,r)vp±{t,r)^0 Vp,qeR' 

/oo 
dx u[ {t, x)up {t, x) = i^"^ S{p-q) 
-oo 



(18) 
(19) 
(20) 
(21) 



/oo 
dx v^g {t, x) Vp {t, x) = i^^ 5{p-q) 
-oo 



(22) 



/oo 
da; u'^g (t, x) Vp (t, x) = Vp, g e R 

-oo 



(23) 



j dr ul± {t, r) Up± (i, x) = = j dr u}^± {t, r) Vp± {t, x) (24) 



y dr vl± (t, r) Vp± {t,x) = = j dr vj^ (t, r) Up± (t, x) (25) 



The electric current density vector is provided by the invariance under U(l) 
phase transformations, according to the Noether theorem, and reads 



cJo(t,r) = g{t,r) = q^^t,x,y)^{t,x,y) (26) 
J,(t,r) = qvF^{t,x,y)-f'^{t,x,y) = - qvp^^ {t,x,y) as^{t, x,y) (27) 
Jy{t,r) = qvF^{t,x,y)-f^^{t,x,y) = qvF^Ht,x,y) a2^{t,x,y) (28) 

where q is the charge. For example, an outgoing particle (massless electron or 
negatively charged neutrino) of momentum p = [p, k) , k ^ , and negative 
charge — e carries the spatially homogeneous current density 

Q = -eul^{t,r)up+{t,r) = ^-j— (29) 

Jx{t;p) = -^exp{-i7rA} 

X (|D,A/2(-^-)|'-iA|Da/2-i(-^-)|') (30) 

X Da/2{-Z-) 1(1-1) VXD^a/2-i{-z+) + c.c. (31) 

Now, since the transverse momentum of a normal mode solution is different 
from zero by definition, it turns out that in the weak field limit A — oo the 
current density components of all the normal modes become exponentially 
small, for any outgoing and/or incoming particle and/or antiparticle normal 
mode massless Dirac spinor. Conversely, for the remaining set of longitudinal 
normal modes with = A; = we immediately find for quasiparticles 

e^E e^E 
g = -T—, Jx = --r-, Jy = 0, VpGR (32) 

nvp n 

while for antiquasiparticles we obviously find 

e^E e'^E 
f?=^, "^- = 1^' VpeR (33) 

This actually means that the only observable nonvanishing current density, 
in the weak field limit, is the longitudinal one, so that the minimum quantum 
conduct ivitj|§ for a graphene sample in the presence of a uniform electric field 



^Notice that in physical units we have e^/h ~ 3.5 x 10'' cm s ^ 4 x 10 ^ 
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is provided by cTmin = 4e^//i ~ 1.4 x 10^ cm ~ 1.5 x 10~^ as indeed 
observed experimentally, which entails Pmax — 6.6 KQ , where the factor 
four is due [6l [7j to the two inequivalent irreducible two-dimensional spinor 
representations, as well as to the so called valley degeneracy. We remark 
that the present value of cxmin is precisely the very same obtained in [1], but 
from a quite different method based upon the Euclidean effective action. 

On the other side, i.e. in the strong field limit A ^ , we readily get for 
graphene quasiparticles 

2e'^E/nhvF , VpeR^Vfc^O . . 

-Ae^E/hvF, \fpeRyk = 

-2e^E/7ih, VpGR^V/c^O , . 

-Ae'^E/h, VpGRVA; = ^ 

Jj^(A = 0) = 0, VpGR^ (36) 

while for antiquasiparticles we obviously find the opposite signs. This means 
that all the quantum states give leading constant contributions, in the 
strong field regime, to the conductivity of a graphene single layer sample, i.e., 
=F4e^/27r/i for the quasiparticle/antiquasiparticle normal modes with k ^ , 
while ^Ae^/h for the corresponding zero modes. This feature nicely explain 
the rapid fall down of graphene resistivity with increasing the charge carriers 
density, viz., their very high mobility. As a matter of fact, if there are N 
electrons in the graphene sample, for instance, then the leading resistivity 
for a sufficiently strong electric field becomes 

-h f 1 forfc = , . 

P ^ 4iV? ^ I 27r for M ^^^^ 

which is vanishing for e.g. a typical density n = 10^^ electrons cm~^ , as 
experimentally observed. 



3 Discussion and conclusions 

In this note we have analysed the exact current density of graphene within 
the pseudo-relativistic effective field theoretic model of a massless Dirac field 
in a 2+1 dimensional space-time under the influence of a uniform background 
electromagnetic field. We have shown that, in the limit of a vanishing electric 
field, a minimum quantum conductivity (Xmin does survive, which is entirely 
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due to the longitudinal zero modes, i.e. with k = 0, and the actual value 
of which is in agreement with the experimental finding. Moreover, it has 
been proved that in the strong field regime all the charged graphene quantum 
states provide a leading, dominant, constant contribution to the conductivity, 
a feature that simply explains the high charge carriers mobility and the rapid 
fall down of resistivity with increasing number of the charge carriers. This is 
the main prediction of the massless, pseudo-relativistic planar QED effective 
model for graphene, in the presence of a constant homogeneous electic field. 
In turn, the related minimum quantum conductivity perfectly reproduces its 
measured values, at variance with many previously obtained different results 
[3j based upon the Kubo formula. The ultimate reason for this might be that 
our exact solutions are truly non-perturbative, viz., they do not reproduce 
the free field spinor when E ^ , whilst the Kubo approach is a perturbative 
linear response approximation. 
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